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ABSTRACT 

We study the evolution of the stellar mass density for the separate families of bulge-dominated and disk- 
dominated galaxies over the redshift range 0.25 < z < 1.15. We derive quantitative morphology for a statisti- 
cally significant galaxy sample of 1645 objects selected from the FORS Deep and the GOODS-S Fields. We 
find that the morphological mix evolves monotonically with time: the higher the redshift, the more disk systems 
dominate the total mass content. At z ~ 1, massive objects (M* > 7 x 10 10 M Q ) host about half of the mass 
contained in objects of similar mass in the local universe. The contribution from early and late type galaxies 
to the mass budget at z ~ 1 is nearly equal. We show that in situ star formation is not sufficient to explain the 
changing mass budget. Moreover we find that the star formation rate per unit stellar mass of massive galaxies 
increases with redshift only for the intermediate and early morphological types, while it stays nearly constant 
for late-type objects. This suggests that merging and/or frequent accretion of small mass objects has a key role 
in the shaping of the Hubble sequence as we observe it now, and also in decreasing the star formation activity 
of the bulge-dominated descendants of massive disk galaxies. 

Subject headings: cosmology: observations — galaxies: morphology — galaxies: formation — galaxies: 
evolution — galaxies: mass function — surveys 



1. INTRODUCTION 

iLillv etafl fi9%). and soon after lMadau et al.l lll996l) . pro- 
duced the first estimates of the cosmic star formation his- 
tory. Since then, the growing number of deep extragalac- 
tic surveys has allowed galaxy evolution to be tackled in 
more and more detail, to higher and highe r redshifts (e.g. 
i Giavalisco et al.ll2004l iGabasch et alJl2004bt iBouwens et al.1 
120041 lluneau et al.ll2005l) . All these studies point towards a 
substantial amount of star formation at early cosmic epochs. 

More recently, a number of stu dies, mainly ba sed on 
NIR s elected surveys lik e 2M ASS llCole et all 120011) . MU- 
NICS jDrorv et all2004al) . K20 (IFontana et al 120041). F IRES 
jRudnick et alJ 120031). HDF JPickinson et alJ 12003) and 
GOODS+FDF (iDrorv et all2005l) . have been able to measure 
directly the stellar mass density up to high redshifts. The two 
independent approaches are in a remarkably good agreement, 
and suggest that about half of the present-day stars was al- 
ready in place at z w 1, when the Universe was half of its 
present age. 

The assembly of the stellar mass through cosmic time is 
a crucial test for models of galaxy formation and evolution 
( Kauffmann & Chariot 1998). Such models aim at linking the 
hierarchical growth of dark matter structures to the observed 
galaxy properties, by means of simplified prescriptions for 
the formation of baryonic systems in d ark matter haloes (e.g. 
iCole et al.l2000tlDe Lucia et al.l200l. 

An effective way to put sensitive constraints on the mod- 

** Based on observations collected at ESO, Chile (ESO Programmes 63. 0- 
0005, 64.O-0149, 64.O-0158, 64.O-0229, 64.P-0150, 65.O-0048, 65.O-0049, 
66.A-0547, 68.A-0013, 69.A-0014 and LP168.A-0485), and on observations 
made with HST/ACS (GO Proposals 9425 and 9502), obtained from the 
ESO/ST-ECF Science Archive Facility. 
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els is to understand where the stars were located, specifically 
what was the morphology of their host galaxies, at differ- 
ent look-back times. In fact, this permits to directly probe 
how galaxies assembled their stars, and how their morphol- 
ogy (i.e., at least with some approximation, their dynamical 
status) evolves. Since merging is driving both mass assem- 
bly and dynamical evolution in a hierarchical scenario, these 
studies offer a direct probe of the role of this process in galaxy 
evolution. 

Because of the lack of sufficient angular resolution, galax- 
ies have been often classified, both at low and high red- 
shift, as early or lat e types based on their broad band colors 
(Baldr v et al1l2004l IFontana et al.ll2004l) . However this kind 
of approach suffers, especially at high redshift, from the ob- 
vious drawback that a disk galaxy populated by an old stellar 
population would be classifie d as an early-type object (and 
vice versa). iBell et all l|2Q03) used the concentration param- 
eter to discriminate between early and late type objects in a 
complete sample extracted from local surveys. They found 
that in the local Universe the transition mass, i.e. the mass at 
which disks become dominant in the relative contribution to 
the total stellar mass function, is rs 5 x 10 10 M Q . 

In this work we rely on a deep, complete and automatically 
morphologically classified sample to study the contribution 
of galaxies of different morph ologies to the reds hift evolu- 
tion of the stellar mass den sity. iBundv et alJ J2005I) . and also 
Brinchmann & Ellis (2000), have carried out a work qualita- 
tively similar to the one we present here, but based on shal- 
lower samples and visual morphological classification. 

This Letter is organized as follows: in ^2] we discuss the 
dataset on which this work is based, in ^3]we present the sur- 
face brightness profiles modeling, and in |4] our results on the 
evolution of mass functions, mass densities and specific star 
formation rates (SSFR, i.e. star formation rate per unit stellar 
mass) split by morphological type. Finally, in g5]we draw our 
conclusions. 

We use AB magnitudes and adopt a cosmology with 
n M =03, n A =0.7, and H =70 kms^Mpc" 1 . 
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FIG. 1 . — Correlation between the Sersic index and the visual classification 
T value for the GOODS-S subsample. Points are mean values. Solid error 
bars are the errors on the means, dashed error bars show the rms values. 



2. GROUND-BASED DATA, PHOTO-Z AND M»/L RATIOS 

This study is based on photometric catalogs deriv ed for 
the FO RS Deep Field (FDF, iHeidt et all2003llGabasch et all 
2004a) and part of the GOODS-South Field (Salvato et al., 
in preparation). The two fields cover approximately the same 
area (39.8 arcmin 2 for FDF and 50.6 arcmin 2 for GOODS-S). 
The FDF photometric catalog covers the UBgRIzJKs broad 
bands, plus an intermediate band centered at 834nm. The I- 
selected catalog used here lists 5557 galaxies down to /; !m = 
26.8. For the GOODS-S field, our K-band selected catalog 
contains 3240 gal axies, and it consists of UB VRIJHK broad- 
band photometry (lArnouts etalJl2001t ISchirmer et"aT]l2003L 
Vandame et al., in preparation). Number counts match the 
literature values down to Ku m « 25. 

Photometric redshift s are d erived using the method de- 
scribed in Bend er et alJ 112001 ). Their accuracy is 3% and 5% 
for FDF and GOODS-S respec tively, with m 1% out liers for 
both fields (for more details see Gabasch et al. 2004a). 

For all galaxies in both catalogs, the M/L ratios were esti- 
mated with a log-likelihood based SED fitting technique, us- 
ing a library of S EDs built with th e lBruzual &Charlotl (120031) 
code, assuming a Salpeter ( 1955) IMF. The procedure used is 
described in detail in Drory et. jd] l|20p4b[). We use here the 
same mass catalog as in Drorv et al. (2005), so we refer the 
reader to this paper for more details. 

3. HST IMAGING AND THE MORPHOLOGICAL ANALYSIS 

Both the FDF and GOODS-S fields were imaged with 

the ACS camera on board the HST. The FDF was imaged 

in the broad-band F814W filter, with 4 WFC pointings 

of 40 minutes exposure each, reaching a 10cr limit of 26 

mag. The data reduction was performed with the standard 

CALACS 5 pipeline, and the combined final mosaic was pro- 

i 1 1 1 I 

duced with the multidrizzle package (Mutchler et al. 2002). 

The GOODS-S field was imaged in four diff erent filters 
(BVIz ) as part of the ACS GOODS legacy JGiavalisco et alJ 
2004); the 10cr limit in the F775W band is 26.5 mag. 

We make use of the publicly available packages GIM2D 
JSimardet all 19991) and GALFIT llPen g et al.l2002l) to fit PSF 
convolved Sersic ( 1968) profiles to the two-dimensional sur- 
face brightness of each object, down to a limit of F814W=24 

5 www.stsci.edu/hst/acs/analysis 



(FDF) and F775W=24.5 (GOODS-S). The PSFs used to con- 
volve the profiles were obtained for each individual tile by 
stacking about 10 high S/N isolated stars. The results from 
the two different codes are in excellent agreement, thus con- 
firming the robustness both of our modeling and of the choice 
of flux limits. 

More details on the reduction and the data quality assess- 
ment, as well as on the whole fitting procedure and relative 
validation, will be described in Pannella et al. (in prepara- 
tion). In order to further check our morphological classifi- 
cation, _vj^ual^orphology\vas_performed by UH according 
to the Ide V aucouleurs et al. (1991) classification scheme. A 
tight correlation is obtained between the average visual and 
automated classifications, parameterized by the morphologi- 
cal type T and the Sersic index n ser , respectively (see Fig.Q. 
Using the visual or automated classification does not signifi- 
cantly affect the results presented in this study. 

We split our sample, according to the Sersic index, in early- 
type (n > 3.5, ss T < -3), intermediate (2 < n < 3.5, « -3 < 
T < 2) and late-type (n < 2, ss T > 2). We evaluated the 
effect of morphological k-correction on this coarse classifica- 
tion scheme by also performing the surface brightness model- 
ing on B and z-band images for a subsample of the GOODS-S 
galaxies. We found that our broad classification is robust for 
objects at z < 1.15. Thus we restrict our analysis to z < 1.15. 

Finally, we estimate the mass completeness of the two 
(FDF and GOODS-S) catalogs as the mass of a maximally 
old stellar population 6 having, at every redshift, an observed 
magnitude equal to the catalog completeness magnitude (e.g. 
iDickinson et alJ2003l) . The final catalog contains 1645 galax- 
ies. 



4. THE EVOLUTION OF THE MORPHOLOGICAL MASS FUNCTION, 
TOTAL MASS DENSITY AND SSFR 

We performed extensive Monte-Carlo simulations to take 
into account the effect of mass uncertainties (« 0.2dex) on our 
results. Ten thousand simulations of the mass catalog were 
generated, perturbing each mass within a gaussian of sigma 
equal to its error. Unless stated differently in the relevant Fig- 
ure captions, we use the median values of the Monte-Carlo 
simulations in all figures. Error bars ta ke into accoun t both 
poissonian errors on the median counts (Gehrels 1986), and 
16-84''' percentile values of each distribution. 

In Fig. |3 we show the VN max corrected stellar mass func- 
tions, split by morphological type, in three redshift bins. The 
redshift intervals were chosen to have comparable numbers 
of objects in each bin. A dotted vertical line shows the mass 
completeness for each redshift bin. We find that, at all red- 
shifts, the late-type objects dominate the low-mass end of 
the mass function, while the early- and intermediate-type ob- 
jects dominate the high-mass end. In addition, Fig.[2]shows 
that the relative contribution of disks to the high-mass end of 
the mass function increases with redshift. This suggests that: 
i) the morphological mix at the high-mass end evolves with 
redshift; ii) the transition mass increase s with redshift. Th e 
latter finding was already suggested in Bundv et al. (2005), 
but its significance was hampered by their mass complete- 
ness limit. In our case, the transition seems to happen well 
above the estimated mass completeness. As Fig. |2| shows, at 
Z ~ 1 the disks' and bulges' contributions become compara- 
ble at w 1 x 10 u Mq. In order to better explore the evolu- 

6 We used a dust-free, passively evolving stellar population model, ignited 
by an instantaneous burst of sub solar metallicity at z = 10. 
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FIG. 2. — The mass function split by morphological types (early in red, intermediate in green, and late in blue) in three redshift bins. Black symbols refer to 
total values. The vertical axis is in units of [hy Mpc~ 3 dex -1 ]. A vertical dotted line indicates the mass completeness limit in each redshift bin. Dashed colored 
lines are intended to guide the eye. Solid black lines show the local mass function determinations from Cole et al. 1 2001 ) and Bell et al. 1 2003). 




FIG. 3. — Top panel: The redshift evolution of the ratio of bulges and 
disks number densities. Error bars show the 5''' and 95' 1 ' percentiles of the 
ratio distributions. Bottom panel: Stellar mass densities for objects with 
M, > 7 X 10 10 Mq (color coding according to galaxy morphology, as in 
Fig.Gj. Points are slightly shifted in redshift for clarity. The median numbers 
of objects in each class are also labeled. Error bars account for poissonian 
eiTors, mass uncertainties, and cosmic variance (for the total values only) as 
described in the text. Dotted and dashed lines indicate the 3<r range esti- 
mate obtai ned at redshift 0, w ith the same mass cut and IMF, from C ole et afl 
l200l and lBell et ai]<20Q3T) . respectively. 



tion of the morphological mix at the high-mass end, we cut 
the three subsamples at the same common mass completeness 
(M* > 7 x 10 I0 M Q ) hereafter. This guarantees a fair com- 
parison between the three redshift bins. The upper panel of 
Fig.|5]shows the behavior of the ratio of the number densities 
of bulges and disks. At increasing redshift, this ratio is found 
to systematically decrease with respect to the no-evolution 
hypothesis at more than 3<r significance. In the bottom panel 
of Fig. |3] we show the stellar mass density split by morpholog- 
ical type at different redshifts. The error estimates of the total 
v alues include the cosmic variance contribution, estimated as 
in S omerville et al.H2004l) . 

In agreement with previous studies (see §1), we find that at 
Z ~ 1 massive objects host almost half of the stellar mass con- 
tained in similarly massive objects at z = 0. In addition, we 



find that the contribution to the total mass budget from early 
and late type galaxies is almost equal at z ~ 1, but strongly 
evolves with redshift: there is a mass pouring from disk sys- 
tems to bulge-dominated objects. Our findings require a sce- 
nario in which massive objects almost double their mass from 
redshift 1 to 0, and become more and more bulge-dominated 
systems as time goes by. 

A possibility is that in situ star formation increases 
the mass and secular evoluti on modifies the morphology 
jKormendv & Kennicutll 120041) . Under simplified assump- 
tions, we can estimate whether the star formation rates (SFR) 
of the galaxies in our three redshift subsamples can account 
for a doubling of the stellar mass by z = 0. Following 
Mada u et ail (119981) . the SFR of individual objects can be 
estimated from the restframe UV luminosity as SFR2800 = 
1.27 x 10~ 28 x L2800 in units of Moyr" 1 , where the constant 
factor is computed for a Salpeter IMF. We apply a dust atten- 
uation correction of A2800 = lmag to the whole sample. This 
med ian value, assum ed not to depend on redshift, is derived 
as in lGabaschl ( 12004) by comparing the total stellar mass den- 
sity and the integral of the SFR density at different look-back 
times. Fig. |4] shows the median SSFR values for the three 
redshift bins and for the different morphological types, cut at 
the same mass completeness as described above. Since it is 
not possible to know the individual star formation histories, 
we derive the median SFR decline with cosmic time from the 
massive galaxies subsample, and we assume that all the ob- 
jects in our sample follow this median star formation history. 
In this way we can estimate, at every redshift, the minimum 
SSFR that enables an object to double its mass by z = 0. These 
minimum SSFRs are shown in the plot as horizontal lines at 
the last two redshift bins considered. It appears that only a 
small fraction of the z ~ 1 sample would be able to double its 
mass. 

An alternative and m ore likely solution (see also 
Brinchmann & Ellis 2000) is a continuous accretion of 
smaller mass galaxies, and possibly the merging of massive 
disks, which can account for both the changing total mass 
budget and the morphological mix evolu tion. 

In the past a number of works (e.g. [ Brinchmann & Ellisl 
l2000tlFeulner et all2005albHBauer et all2005l) have focused 
on the SSFR redshift evolution and consistently found an in- 
crease of the SSFRs with redshift. Hence the bulk of star for- 
mation in massive galaxies was pushed to z > 2, in agreement 
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FIG. 4. — The median SSFR for different morphological types. The 
coloured symbols (same colour coding as in previous figures) show the me- 
dian values in both stellar mass and SSFR. Dashed error bars represent the 
median 10'* and 90''' percentiles of the Monte-Carlo realizations, while the 
solid error bars show lcr uncertainties on the median values. Tilted lines are 
constant SFR values (from the bottom 1,5,40 Mqjt -1 ) while the horizontal 
lines mark the doubling mass lines (only for the last two redshift bins - see 
text for details). 



with the downsizing scenario (ICowie et al 

Here for the first time we explore the SSFR evolution split 
by different morphological types. We find, for objects with 
M,>7x 10 10 M Q , a different behavior in the SSFR redshift 
evolution for the different morphological types. The median 
SSFR of late-type objects shows only a mild, if any, evolution 
up to redshift 1, while the median values of early and inter- 
mediate types are clearly increasing with redshift. Since these 
latter populations are dominanting in numbers, they drive the 
redshift evolution of the whole massive galaxies distribution. 
This finding, together with the previous arguments, suggests 
that the morphological evolution of massive disk galaxies to- 
ward more and more bulge dominated systems is accompa- 
nied by a decrease in the specific star formation rate of their 
e nd products. 

Bell et al. (2005) show that a large part of the star formation 
in massive galaxies at z ~ 0.7 is due to normal spiral galax- 
ies. Moreover, they conclude that turning off star formation in 



these objects is responsible for the observed decline of cosmic 
star formation therafter. We find a good agreement at z ~ 0.7 
with their result but also a hint that the decline of star for- 
mation in massive galaxies is linked more to the evolution of 
bulge-dominated galaxies. 

5. CONCLUSIONS 

In this Letter we have studied the contribution of disks and 
bulges to the evoluti on of the s tellar mass function up to z ~ 1 . 

We agree with Bundv et al. (2005) that the transition mass, 
i.e. the mass at which disks become dominant in the relative 
contribution to the total mass function, increases with redshift. 
We estimate that at z ~ 1 the transition mass is up to a factor 
2 larger than its measured local value. 

We show that the morphological mix evolves with redshift. 
At z ~ 1 early and late type galaxies contribute nearly equally 
to the total mass budget in massive objects (i.e. with M* > 
7 x 10 10 Mq). We suggest that merging events must play a 
key role in the mass pouring from disks to bulges. 

We find a different behavior of the SSFR, i.e. the star for- 
mation rate per unit stellar mass, as a function of redshift for 
the different morphological types. The median SSFR of late- 
type objects shows almost no evolution up to z ~ 1. Con- 
versely, median SSFRs for early and intermediate types in- 
crease systematically. Since these latter morphological types 
are dominating in numbers, they drive the total SFR evolution 
of massive galaxies. This suggests a scenario where the mor- 
phological evolution of massive disk galaxies through merg- 
ing is followed by the decrease of the star formation in their 
bulge-dominated descendants, maybe after a burst of star for- 
mation that exhausts the available gas. 

We wish to thank the referee for valuable comments that 
pushed us to strengthen our results and to improve their pre- 
sentation. M.R is glad to thank G. De Lucia, D. Pierini, 
G. Rudnick and V. Strazzullo for a careful reading of the 
manuscript and for many insightful comments and sugges- 
tions. This work was supported by the Deutsche Forschungs- 
gemeinschaft through the SFB 375 grant. 
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